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We calculate the equation of state of dense hydrogen within the chemical picture. Fluid variational theory is 
generalized for a multi-component system of molecules, atoms, electrons, and protons. Chemical equilibrium 
is supposed for the reactions dissociation and ionization. We identify the region of thermodynamic instability 
which is related to the plasma phase transition. The reflectivity is calculated along the Hugoniot curve and 
compared with experimental results. The equation-of-state data is used to calculate the pressure and temperature 
profiles for the interior of Jupiter. 



1 Introduction 

The equation of state (EOS) of hydrogen and helium at high pressures is of great relevance for models of the 
interior of giant planets and other astrophysical objects as well as for inertial confinement fusion experiments. 
For detailed calculations accurate knowledge of the EOS over a wide range of densities and temperatures is 
needed. Especially, in the range of warm dense matter with high densities characteristic for condensed matter 
and at temperatures of a few eV the EOS is crucial for modelling giant planets. This region is challenging 
for many-particle theory because strong correlations dominate the physical behavior. Progress in shock-wave 
experimental technique has allowed to study this region only recently. 

To probe the EOS, experimental investigations were performed statically with diamond anvil cells or dynami- 
cally by using shock waves, see fl] for a recent review. The experimental data indicate that a nonmetal-to-metal 
transition occurs at about 1 Mbar which is identified by a strong increase of the conductivity [2] and reflectiv- 
ity Q. Some theoretical models yield a thermodynamic instability in this transition region, the plasma phase 
transition (PPT) [4, 5 6 7, 8 |, which would strongly affect models for planetary interiors and the evolution of 
giant planets ll9l [T0l[TT1l . After a long period of controversial discussions, new results of shock wave experiments 
on deuterium support the existence of such a PPT fl2l . This fundamental problem of high-pressure physics will 
also be studied with the FAIR facility at GSI Darmstadt within the LAPLAS project, see ifHl [1411 . 

In this paper we present new results for the EOS of dense hydrogen within the chemical picture. We treat the 
reactions pressure dissociation and ionization self-consistently via respective mass action laws. We identify the 
region of thermodynamic instability and calculate the phase diagram as well as the reflectivity in order to verify 
the corresponding nonmetal-to-metal transition. The EOS data is used to model the interior of Jupiter within a 
three-layer model. The agreement with astrophysical constraints such as the core mass and the fraction of heavier 
elements can serve as an additional test of the theoretical EOS. 

2 Equation of state for dense hydrogen 

Warm dense hydrogen is considered as a partially ionized plasma in the chemical picture. A mixture of a neutral 
component (atoms and molecules) and a plasma component (electrons and protons) is in chemical equilibrium 
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with respect to dissociation and ionization. The EOS is derived from an expression for the free energy of the 
neutral (i*b) and charged particles (F±), see lfT31[T6l : 

F(T,V,N) = F +F±+F pol . (1) 

The first two terms consist of ideal and interaction contributions and can be written as F = F^ d + F in * and F± = 
F± + F± lt . Fp i contains interaction terms between charged and neutral components caused by polarization ifTTl . 

Applying fluid variational theory (FVT), the EOS is determined by calculating the free energy Fq"*(T, V, N) 
via the Gibbs-Bogolyubov inequality |[T8l . This method has been generalized to two-component systems with a 
reaction |fl9l l20l |2D so that also molecular systems at high pressure can be treated where pressure dissociation 
occurs, e.g. H2 ^ 2H for hydrogen. In chemical equilibrium, /xh 2 = 2/ijj is fulfilled, and the number of 
atoms and molecules can be determinded self-consistently via the chemical potentials /i c = (dF/dN c )T- The 
effective interactions between the neutral species are modeled by exp-6 potentials, and the free energy of a multi- 
component reference system of hard spheres has to be known; for details, see |[T9ll20ll22l . 

The charged component is treated by using efficient Pade approximations for the free energy developed by 
Chabrier and Potekhin l23ll . The coupling with the neutral component occurs via the ionization equilibrium, 
H^e+p. In chemical equilibrium, the relation /zh = Me + Mp determines the degree of ionization. 

Since atoms and molecules are particles of finite size there is an additional interaction between the charged 
component and the neutral fluid. According to the concept of reduced volume, point-like particles cannot pene- 
trate into the volume occupied by atoms and molecules. This leads to a correction in the description of the ideal 
gas of the charged component 112411251 so that the ideal free energy of protons and electrons F± is dependent on 
the reduced volume V* = V ■ (1 — rj), 

F±(T, V*,N) = N ± k B T ■ (2) 

where r\ is the ratio of the volume which cannot be penetrated by point-like particles to the total volume. It is 
derived from hard sphere diameters obtained within the FVT self-consistently. The free energy density f±'* is 
given by Fermi integrals which take into account quantum effects. In order to avoid an intersection of pressure 
isotherms, which is important for modelling planetary interiors, a minimum diameter d m ; n has been introduced. 
It was determined starting at low temperatures where it remains almost constant up to 15.000 K, then it increases 
up to 20.000 K and remains constant again for higher temperatures, see Fig.Q] These values are in the range of 
the results for the diameter of the hydrogen atom derived from the confined atom model [26|. 
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Fig. 1 Minimum diameter for expanded particles (atoms, molecules) introduced within the reduced volume concept. 

Consequently, the reduced volume concept changes the chemical potential of each component drastically at 
higher densities and results in pressure ionization. This is due to the fact that additional terms appear in the 
chemical potential, which is the particle number derivative of the free energy, and thermodynamic functions of 
degenerate plasmas are very sensitive to changes in density. 

This current model FVT + includes all interaction contributions to the chemical potentials, thus being a gener- 
alization of earlier work 11221 where only ideal plasma contributions have been treated (FVTj^). 
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Fig. 2 Composition of dense hydrogen for 20.000 K (left) and 50.000 K (right). 



In Fig. |2] the composition of hydrogen derived from the present approach is shown for two temperatures. 
Hydrogen is an atomic gas at low temperatures (left) and low densities. With increasing densities molecules 
are formed due to the mass action law. Pressure dissociation and ionization can be observed in the high-density 
region. The nonideality corrections to the free energy force a transition from a molecular fluid to a fully ionized 
plasma. At higher temperatures (right) the formation of molecules is suppressed and pressure ionization becomes 
the dominating process. At low densities and high temperatures a fully ionized plasma is produced due to thermal 
ionization. 

We show pressure isotherms over a wide range of temperatures and densities in Fig. [3] At low densities the 
system behaves like a neutral fluid. Between densities of 10~ 3 g/cm 3 and 10 -1 g/cm 3 nonideality corrections to 
the free energy of atoms and molecules lead to a nonlinear behavior of the isotherms. For still higher densities a 
phase transition occurs which is treated by a Maxwell construction. The thermodynamic instability vanishes with 
increasing temperatures, and the critical point is located at 16.800 K, 0.35 g/cm 3 , and 45 GPa. 




Fig. 3 Pressure isotherms for dense hydrogen. 



The critical point and the related coexistence line are shown in Fig.|4]and compared with results of other EOS. 
The critical point itself lies within the range of other predictions, whereas the coexistence line is lower than most 
of the other results. For a comparison of data concerning the PPT, see TableQ] 

New shock-wave experiments [12| imply that a PPT occurs in deuterium at densities of 1.5 g/cm 3 and a 
coexistence pressure of about 1 megabar. Each of these values is twice as high as evaluated in the recent model. 
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Fig. 4 Phase diagram for dense hydrogen. Present results of the FVT + (red) are compared with other predictions for the 
PPT: SC El[5l, RK (27), MH HE), ER |29|, SBT |30], RRN ED, BEF |32|, MCPB (33). 
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Table 1 Theoretical results for the critical point of the hypothetical plasma phase transition (PPT) in hydrogen which was 
predicted by Zeldovich and Landau | 37] and Norman and Starostin (38). 



3 Conductivity and reflectivity 

The PPT is an instability driven by the nonmetal-to-metal transition (pressure ionization). We calculate the 
electrical conductivity as well as the reflectivity by applying the COMPTRA04 program package [39] [40) in 
order to locate this transition in the density-temperature plane. 

Optical properties are calculated within the Drude model. The reflectivity R(oj) is given in the long-wavelength 
limit via the dielectric function e(uj) which is determined by a dynamic collision frequency v(u) or, alternatively, 
by the dynamic conductivity <j(oj) ATI : 



a{uj) = o-(0) 



y/e$*)-l 

2 

LU 

1 - 



l 

r m' = 1 + 

U) [OJ + 1V(UJ)\ EqUJ 

-, -1 

— a(0) 



1 - 



(3) 
(4) 
(5) 



yj n e e 2 / (eoTO e ) is the plasma frequency of the electrons. 



The reflectivity was determined along the Hugoniot curve and is compared with experimental results [ 3 1 and 
those of the earlier model FVTj^ fl42ll in Fig. [5] The results of the current model show a much better agreement 
with the experiment. The characteristic and abrupt rise with increasing pressure was reproduced more accurately. 
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This drastic increase appears due to pressure ionization in the vicinity of the criotical point of the PPT. As a 
result, the reflectivity advances from very low values to metallic-like ones almost instantly. 
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Fig. 5 Reflectivity of dense hydrogen within the models FVT + and FVTj^ along the Hugoniot curve in comparison with 
experiments (3). 



4 Planetary interiors 

Modelling the interiors of giant planets and comparison with their observational parameters offers an alternative 
tool besides laboratory experiments of probing the EOS of the components the planets are predominantly made 
of. Giant planets such as Jupiter and Saturn consist mainly of hydrogen and, in decreasing order, of helium, 
water and rocks, covering a wide range of pressures and temperatures. Independently from the H-EOS used for 
modelling, the simplest interior structure that is compatible with the observational constraints requires at least 
three homogenous layers with a transition from a cold molecular fluid in the outer envelope to a pressure ionized 
plasma in the deep interior and a dense solid core of ices and rocks. A solid core may be explained as a result 
of the formation process and the seperation into two fluid envelopes with different particle abundances by an 
existence of a PPT as provided by the FVT+ EOS. The constraining observational parameters are the total mass 
of the planet M, its equatorial radius R eq , the temperature T at the outer boundary, the average helium content 
Y , the period of rotation u> and the gravitational moments J 2 , J4 , Jq . From measurements of the luminosity it 
has been argued |43| that the temperature profile should be adiabatic. For a given EOS, the interior profiles of 
pressure P and density p are calculated by integration of the equation of hydrostatic equilibrium 

jjLj V ? P(r, 4 . V, (o I <?r'0- { + \ U V sto P) (6) 

along an isentrope defined by the outer boundary. The first term on the right hand side of eq. ©is the gravitational 
potential and the second term the centrifugal potential assuming axialsymmetric rotation. We apply the theory 
of figures [44 1 up to third order to solve this equation and to calculate the gravitational moments. They are 
defined as the coefficients of the expansion of the gravitational potential into Legendre polynomials, taken at the 
outer boundary. Being integrals of the density distribution weighted by some power of the radius, they are very 
sensitive with respect to the amount and distribution of helium and heavier elements within the planet. 

In accordance with previous calculations B31 1461 . mixtures of hydrogen with helium and heavier elements 
have been derived from the EOS of the pure materials via the additive volume rule. It states that the entropy of 
mixing can be neglected. 
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Assuming a three-layer structure, we present results for Jupiter for the profiles of temperature, density, and 
pressure along the radius in Fig. fusing two different H-EOS, the standard Sesame table 5251 for hydrogen [47] 
and the FVT + model presented above. 

profiles inside Jupiter 
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Fig. 6 Profiles of temperature, density, pressure along the radius within Jupiter using two different H-EOS, FVT + (solid) 
and Sesame 5251 (dashed). 

The profiles of temperature appear very similar, meaning a small uncertainty about the real profiles. Contrary, 
the density and pressure profiles exhibit more differences and require some explanation. In the fluid part of 
Jupiter, the presence of a PPT leads to a jump in density between the envelopes. Since the gravitational moments 
as integrals over the density have to be the same for both H-EOS, the density profile of a H-EOS with PPT has to 
be smaller in the outer envelope and larger in the inner envelope. The different size and composition of the core 
for these specific H-EOS are a consequence of their different compressibility in the regime of pressure ionization 
at about 1 Mbar, where the gravitational moments are most sensitive to the density distribution. 

In case of a stiff H-EOS like Sesame, a larger amount of heavy elements is needed in the two fluid envelopes 
to compensate for the smaller hydrogen density at a given pressure. As a result, this material is added to the 
well-known density-pressure relation of degenerate electrons in the deep interior, leaving less material for the 
core. Thus, in case of the Sesame-EOS, the amount of heavy elements becomes with 10% very large and an 
unlikely solution with a very small core of light material (e.g. water) can be found. 

In case of the FVT+ EOS which is more compressible than the Sesame EOS at about 1 Mbar, the helium 
content is below the value of 27.5% for the protosolar cloud in order to reproduce the lowest gravitational moment 
J2. Furthermore, the next gravitational moment J4 cannot be reproduced correctly because the transition to the 
metallic envelope occurs already at about 90% of the radius and, thus, at too low densities. For opposite reasons, 
both the Sesame and FVT+ EOS applied in a three-layer model of Jupiter are not compatible with all of the 
observational constraints. While Sesame is probably too stiff, the FVT + model is likely too soft in the WDM 
region at about 1 Mbar. 

5 Conclusions 

In this paper, we have extended the earlier chemical model FVTi to calculate the EOS of dense hydrogen. 
The current model FVT + includes nonideality corrections to the free energy of each commponent of the partially 
ionized plasma. We have shown results for the composition and the thermodynamic properties of dense hydrogen. 
The PPT was located in the phase diagram, its critical point coincides with earlier results. Furthermore, we 
have determined optical properties such as reflectivity and conductivity, within linear response theory using the 
program package COMPTRA04. The calculated reflectivity along the experimental Hugoniot curve shows a 
good agreement with the experiments. However, application of the FVT+ EOS to the interior structure of Jupiter 
indicates that the behavior at about 1 Mbar is probably too soft. The same conclusion can be drawn from a 
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comparison with shock-wave experiments that indicate the existence of a PPT lfl2l . FVT + predicts the PPT at 
too low pressures as well as at too low densities. Further efforts to solve this problem, especially concerning the 
reduced volume concept, are necessary. 
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